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ABSTRACT 

We consider a coplanar system comprised of a massive central body (a star), a less 
massive secondary (a planet) on a circular orbit, and a test particle on a bound orbit 
exterior to that of the secondary. The gravitational pull exerted on the test particle 
by the secondary acts as a small perturbation, wherefore the trajectory of the particle 
can be described as an ellipse of a precessing perihelion. While the apsidal motion 
is defined overwhelmingly by the Newtonian portion of the secondary's gravity, the 
post-Newtonian portion, too, brings its tiny input. We explore whether this input may 
be of any astrophysical relevance in the next few decades. We demonstrate that the 
overall post-Newtonian input of the secondary's gravity can be split into two parts. 
One can be expressed via the orbital angular momentum of the secondary, another via 
its orbital radius. Despite some moderately large numerical factors showing up in the 
expressions for these two parts, the resulting post-Newtonian contributions from the 
secondary's gravity into the apsidal motion of the test particle turn out to be small 
enough to be neglected in the near-future measurements. 

Key words: celestial mechanics - gravitation - planets and satellites: general - 
methods: analytical 



1 INTRODUCTION 

The perihelion shift of Mercury gave an experimental ev- 
iden ce for the t heory of general relativity a century ago, 
e.g. (IWiU 1 1200611 . This relativistic effect is caused by the 
solar gravity, which is 10~* smaller than the Newtonian 
ones. Substantial improvements including C assini mission 
have been brought into the ephemerides, e.g. (|Pitieva II2008I . 
I2OI0I ). where the post-Newtonian effects in the solar system 
have been considered n umeric ally in a co mbined manner, 
e.g. (|Maindl and Dvora k 1994: Simon et a l. 1994). 

We should note that the post-Newtonian effects of plan- 
etary gravity on a light-like trajectory are of astrophysical 
relevance. Two decades ago, the bendi ng angle of light by a 
plane tary mass was observed by VLBI (|Treuhaft and Lowe I 
I1991I ). It is interesting to investigate the post-Newtonian ef- 
fects of planetary gravity on a particle motion. 

We may insist, before doing explicit calculations, that 
even the post-Newtonian effects of the jovian gravity should 
be 10~^ smaller than the solar ones and thus they should be 
negligible. This argument seems reasonable, unless a large 
numerical factor comes in. For instance, 10~^ multiplied by 
say 67r^ is nearly five percents. It is not the best idea to 
use numerical simulations for the purpose of discussing the 
smaller post-Newtonian effects of planetary gravity sepa- 
rately from other major effects. 



The main purpose of this paper is to clarify whether 
a large numerical factor is involved in the post-Newtonian 
effects of planetary gravity on the perihelion shift. For this 
purpose, we propose a model calculation that is easily ac- 
cessible. 

This paper is organized as follows. In section II, analyti- 
cal calculations are carried out to study the post-Newtonian 
effects of planetary gravity on the perihelion shift. The ef- 
fects are evaluated for the solar planets in section III. Section 
IV is devoted to the conclusion. We provide some calcula- 
tions in Appendices. 

Latin indices run from 1 to 3 and Greek ones run from 
to 3. Einstein summation convention is used for both of 
the indices. We take the units of G = c = 1. 



2 POST-NEWTONIAN EFFECTS OF 
PLANETARY GRAVITY 

2.1 Effective metric by an orbiting planet 

In the post-Newtonian approximation, the line element 
for A'^ masses is expressed in an inertial frame as 
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l|Landau and Lifshitz II 19621 : iMisner et al. 1119731 ) 
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where {x'^) — {t,r) for /i = 0,1,2,3, the position of each 
mass niA is denoted as Ra, the relative vectors and distances 
are defined as ta = t - -Ra, Rab = Ra — Rb, rA = [vaI, 
Rab = I -Rab I and r = |r| and the velocity is va- 

The terms in front of dtdx-' in the R.H.S. of Eq. ^ make 
a contribution as shown in the Appendix A. It is caused by 
the orbital angular momentum of both the central object 
and the secondary one. Under the above approximations, 
however, the orbital angular momentum of the central object 
is much smaller than that of the secondary one, so that it can 
be neglected (See Appendix B). This effect is distinct from 
the Le nse-Thirring effect that com es from a single spinning 
object (|Thirring and Lense Ill9ld ). 

We consider three masses rriA {A — 1, 2, 3). See Figure 
[T]for the configuration of the system and our notation. For 
simplicity, approximations are made for the system: 

(i) mi > m2 » ma, 

(ii) £ <C a {R2 ^ r at the third body position), 

(iii) on the same orbital plane (coplanar), 

(iv) a circular motion of the secondary mass, 

where the orbital radius of the secondary object is denoted 
as i and the semi-major axis of the third mass orbit is de- 
noted as a. The approximation (i) implies that the third 
body can be treated as a test particle in the spacetime pro- 
duced by the two-body system of mi and m2. The approxi- 
mation (ii) suggests that the orbital period of the third body 
is much longer than that of the secondary one according to 
Kepler's third law. The approximation (iv) means that the 
primary mass is also moving on circular orbit around the 
common center of mass, because ms is a test mass. 

The post-Newtonian center of mass may affect the mul- 
tipole expansion. Fortunately, the correction by the post- 
Newtonian center of ma ss vanishes for the circular motion 
l|Lincoln and Will Ill99(]| '). 

The slowly-moving body feels the gravitational field 
that is averaged over fast-changing fluctuations caused by 
the moving primary and secondary objects. For discussing 
the motion of the third body, therefore, it is convenient to 
use the time averaging < >t= J dt for the orbital pe- 
riod P. For a circular motion, this averaging is the same as 
the angular one < >^= (27r)~^ J d'd. This simplifles the 
following calculations. The averaging is denoted simply as 
< >. 

Under the approximations (i)-(iv) with the averaging. 



the averaged metric becomes 
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The calculations for deriving this are given in the Appendix 
C, where our estimate shows that the post-Newtonian effects 
stemming from the circular motion of the primary can be 
neglected. 

Let us introduce the polar coordinates as dr'^ + r'^ {d6^ + 
sin^ Bdip^) = Sijdx^dx-' . For a radial coordinate R to be the 
circumference radius, we introduce the isotropic coordinates 
as 
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By using this transformation, we reach the effective 
metric acting on the third mass as 
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where we have used 9 = tt/2 = const, for the coplanar con- 
figuration. Here, we defined the effective Schwarzschild ra- 
dius as rs ~ 2mtot and the effective moment induced by the 
secondary body as 



mar 



(5) 



It follows that the line element given by Eq. (|4]) agrees with 
the Schwarzschild metric on the equatorial plane if Q = 0. 



2.2 Motion in the effective potential 

Equation Q gives the Lagrangian for the third body as 

L = - 



R 



i?3 i?4 



(6) 



where the dot denotes the derivative with respect to the 
proper time of the test mass as d/dr. 

We get constants of motion for this Lagrangian as 



R 



i?3 i?4 



(7) 

J = iiV- (8) 
Substituting these constants into Eq. (|6]) leads to the energy 
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integral as 
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(9) 



Note that this includes both the Newtonian effect with a 
quadrupole potential (expressed as —Q/R^ at the third term 
of Eq. ([9])) and the Schwarzschild case at the first post- 
Newtonian order. 

It is convenient to introduce 
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(10) 



Using dR/dr = {dR/dip){dip /dr) and Eq. (O for Eq. we 
formally obtain 

%)'-Fiu, (U, 

where F{u) is a function of u with coefficients including 
and Q. 
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Figure 1. Schematic figure of the three masses mi, 1712 and m^. 
The mean orbital radius of the secondary mass is denoted as 
The semi-major axis and the eccentricity of the Keplerian orbit 
for the third body are denoted as a and e. 



2.3 Relativistic perihelion shift at 0{Q) 

Let Ta and rp denote the aphelion and perihelion in the 
Keplerian orbit of the third mass (the semi-major axis a 
and the eccentricity e). We define Ua = r^^ and Up = ^, 
which are related with the orbital elements in the Newtonian 
limit as 

1 



a{l + e) 
1 



a{l 



(12) 
(13) 



The perihelion shift for (r = — > Va 
as (|Weinberg Ill972h A(p — 2 JJ^" dip — 2n, where the factor 
2 comes from the fact that the averaged metric by Eq. Q 
is axisymmetric. 

By integrating Eq. (|ll|l , the perihelion shift of the third 
body {ua Up ^ Ua) is expressed at 0{Q) as 

A^Q= (^24+^e" + 0(e*)^7r^. (14) 

Note that r 2> Ra leads to e <C 1 and hence the above 
expression is valid for small eccentricity e <^ 1. Eq. (|14p 
yields the perihelion shift rate as 



24+lIle^ + 0(e- 



Q_7r 
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— e + 0(e ) I — —n, 



(15) 



where n = 2ttP~^ is the mean motion of the third mass. 
The post-Newtonian correction to P could cause the 2PN 
correction to rb. 



Table 1. ro^ and -cuq by Jupiter and its inner planets on the 
Saturn's perihelion, where Eqs. JTSj and IIA3II are used. 



planet 


m2 [kg] 


e [AU] 


■uu]^ [mas/cy] 


■djQ [mas/cy] 


Mercury 


3.30 X 10^^ 


0.387 


-6.09 X 10-^ 


7.50 X 10-8 


Venus 


4.87 X 10^* 


0.723 


-1.23 X 10~^ 


3.86 X lO-'^ 


Earth 


5.98 X 10^* 


1.00 


-1.77 X 10-5 


9.06 X 10-'^ 


Mars 


6.42 X 10^^ 


1.52 


-2.35 X lO"*^ 


2.26 X lO"'^ 


Jupiter 


1.90 X 10^^ 


5.20 


-1.28 X 10-2 


7.80 X 10-3 



momentum among the solar planets. First, we choose the 
secondary object as the Jupiter and estimate the magnitude 
of the jovian effects on the Saturn. Table [1] shows the plane- 
tary effects due to the orbital angular momentum (wl) and 
the induced quadruple moment (ri^g). The orbital angular 
momentum makes a slightly larger effect than the induced 
quadruple moment. 

Next, let us consider inner planets as another example. 
The post-Newtonian effects of the Mercury on the Venus 
are —1.4 x IQ-^ and 8.2 x IQ-* mas/cy, respectively. The 
estimated values are smaller by two (or mo re) digits than 
the uncerta inties in the current ephemerides (jPitieva 1 l2008l . 
I2010l . l201ll ). 

The above corrections to the perihelion shift for thirty 
years amount to ~ 1.3 x IQ-^ [mas/cy] x 30[yr] ~ 4 x 
10-3[mas] for the Saturn and ~ 1.4xl0-3[mas/cy] x30[yr] ~ 
4 x 10^* [mas] for the Venus. Therefore, the post-Newtonian 
effects of the planets on the orbital motion are small enough 
to be negligible even for near future measurements. 



3 SOLAR PLANETS 

Eqs. (|14|l and p5|) involve not huge but moderately large 
numbers such as 247r in the coefficients. Therefore, we make 
an estimation of the corrections for the solar planets. 

The Jupiter has the largest mass and orbital angular 



4 CONCLUSION 

The post-Newtonian effects of planetary gravity do influence 
the perihelion shift. Analytical calculations were carried out 
to show that the post-Newtonian effects of planetary grav- 
ity originate from the orbital angular momentum and the 
orbital radius of a planet. Thereby, we argued that the post- 
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Newtonian effects are small enough to be negligible in the 
solar system even for near future measurements, though a 
moderately large numerical factor is involved. A study be- 
yond the present approximations is left as future work. 

Finally, let us mention a gauge issue of the present re- 
sult. The metric Eq. ([1} is based on the standard PN co- 
ordinates. Namely, the gauge freedom is fixed. The gauge 
choice at the post-Newtonian order makes a d ifference of 
10~^ at most in the s olar system dynamics, e.g. (|Brumberg I 
I1991I : lKlioner1l2003l ). Therefore, we should note that our 
result (e.g. the numerical coefficients) is dependent on the 
coordinates that we used. 
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APPENDIX A: EFFECTS BY THE 
PLANETARY ORBITAL ANGULAR 
MOMENTUM 

Let hoj denote the coefficient of dtdx-' in ds^ expressed by 
Eq. ([1]). The field hoj is generated by massive bodies 1 and 
2. In the weak field approximation, / ip, is known to take the 
form as l|Landau and Lifshitz |[l963 ) 



■oj 



-(n X L) 
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where L denotes the vector for the total orbital angular 
momentum of the central two objects . The perihelion shift 
rate by this metric component reads (jLandau and Lifshitz I 
Il962l ) 
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where L is the spin angular momentum. For our case (mi 2> 
1712), L corresponds to the orbital angular momentum as 
17121^ {mtot / i^y^'^ , where the orbital angular momentum of 
the primary object is negligible for mi ^ m2. Substitution 
of this into Eq. HA2|I leads to the perihelion shift rate as 



4m2^ 
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APPENDIX B: SOLAR AND PLANERARY 
ORBITAL ANGULAR MOMENTA 

The solar and planetary orbital angular momenta [kg ■ 
m s~ ] with respect to the solar system's barycenter are 
1.8 X 10*" (Sun), 9.2 x 10^** (Mercury), 1.8 x 10*° (Venus), 

2.7 X 10*° (Earth), 3.5 x 10^^ (Mars), 1.9 x 10*^ (Jupiter), 

7.8 x 10*^ (Saturn), 1.7 x 10*^ (Uranus) and 2.5 x 10*^ (Nep- 
tune), where the solar orbital angular momentum is caused 
predominantly by outer planets such as the Jupiter and Sat- 
urn. The Sun's orbital angular momentum is much smaller 
than that of the Jupiter (and the Saturn). 



APPENDIX C: AVERAGED METRIC 

Let us explain how to obtain Eq. ((2]). The potentials in the 
metric expressed by Eq. ([T]) are expanded, because r 3> Ra 
at the third-body position. For instance, 

1 _ 1 
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(CI) 
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where we define unit vectors as n = r 
One can use m\R\n\ — — m2-R2«2, because the origin of the 
coordinates is chosen as the center of mass. One can expand 
the Newtonian potential as 

mA _ ratot 
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^ 2 
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Ra 
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where the third body is considered a test mass and hence we 
define the total mass as ratot = mi -I- m2 and cos i9 = n ■ n2 
defines the angle from the direction to the secondary body 
to the third one at the coordinates origin. It is clear that 
the second term in the R.H.S. of Eq. (|C2|I is the quadrupo- 
lar part. We may think that the post-Newtonian center of 
mass may affect the above expansion. Fortunately, the cor- 
rection by the post- Newtonian center of ma ss vanishes for 
the circular motion (|Lincoln and Wiinil99Gl ). 

After being averaged, the Newton-type potential be- 
comes 
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The assumption (i) leads to 
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where £ = i?2i. Hence, we obtain miRf + m2-R2 = rn,2£'^ at 
the lowest order, so that Eq. (|C3|) can be rewritten as 
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Similarly, under the assumptions (i) and 
velocity-dependent part is averaged as 
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where for computing this post-Newtonian part it is sufficient 
to make substitutions as vl — RiD."^ — )• and V2 = i?!^^^ ^ 
. Here, the secondary's angular velocity is denoted 



mtott. 

as fi. Finally, the second-order part in masses is expanded 
as 
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which is not affected by the averaging. 

By combining the above results, therefore, the averaged 
metric is expressed by Eq. ([2]). 
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